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Abstract 

Single nucleotide-polymorphisms (SNPs) are a source of diversity among liuman population, which may be responsible for 
the different individual susceptibility to diseases and/or response to drugs, among other phenotypic traits. Several low 
penetrance susceptibility genes associated with malignant melanoma (IVIIVI) have been described, including genes related to 
pigmentation, DNA damage repair and oxidative stress pathways. In the present work, we conducted a candidate gene 
association study based on proteins and genes whose expression we had detected altered in melanoma cell lines as 
compared to normal melanocytes. The result was the selection of 88 loci and 384 SNPs, of which 314 fulfilled our quality 
criteria for a case-control association study. The SNP rs6854854 in ANXA5 was statistically significant after conservative 
Bonferroni correction when 464 melanoma patients and 400 controls were analyzed in a discovery Phase I. However, this 
finding could not be replicated in the validation phase, perhaps because the minor allele frequency of SNP rs6854854 varies 
depending on the geographical region considered. Additionally, a second SNP (rs6431588) located on ILKAP was found to 
be associated with melanoma after considering a combined set of 1,883 MIVl cases and 1,358 disease-free controls. The OR 
was 1.29 (95% CI 1.12-1.48; p-va/ue = 4xl0"''). Both SNPs, rs6854854 in ANXA5 and rs6431588 in ILKAP, show population 
structure, which, assuming that the Spanish population is not significantly structured, suggests a role of these loci on a 
specific genetic adaptation to different environmental conditions. Furthermore, the biological relevance of these genes in 
MM is supported by in vitro experiments, which show a decrease in the transcription levels of ANXA5 and ILKAP in 
melanoma cells compared to normal melanocytes. 
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Introduction 

Malignant melanoma (MM) is a .stepwise tumor process in 
which normal melanocytes in the basal layer of the epidermis 
acquire genetic aberrations that drive progression to melanoma. In 
recent years, there has been a constant increase in the incidence of 
MM among the world population, which is a worrying fact 
because MM is a highly aggressive, potentially lethal form of 
cancer [1]. Summaries of national melanoma notifications 
provided to the International Agency for Research on Cancer 
(lARC) (2002) demonstrate that the highest reported national 
incidence rates for melanoma occurred in the populations of 
Australia (39:100,000 per year) and New Zealand (34:100,000 per 



year). The next highest national melanoma rates were observed in 
the USA (17:100,000 per year) followed by the countries in 
northern (Denmark, Norway and Sweden) and western European 
countries (the Netherlands and the United Kingdom), with 
incidences of 9-15: 100,000 per year [2,3]. According to data 
published by the Spanish National Epidemiology Carlos III 
Health Institute, the annual incidence of melanoma in Spain is 
close to the adjusted annual incidence for the European 
population of 6.14 per 100,000 population in men and 7.26 per 
100,000 population in women [3]. The predominandy non- 
Caucasian populations of Africa and Asia reported melanoma 
rates less than 3:100,000 per year [2]. Likewise, MM tends to 
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occur more often in people wdth light skin, hair and eyes, who 
seem to be more sensitive to the sun's ultraviolet radiation, 
although the development of this pathology may happen in any 
geographical human group [4]. 

In this regard, epidemiological and genetic studies show that 
both the development of the disease and its evolution are 
determined by individual-specific genetic factors, genetic and 
epigenetic aberrations acquired by the tumor, and by environ- 
mental conditions [5,6,7]. Two high-risk melanoma susceptibility 
genes with large effects but low frequency in the population have 
been described. The best-established high-risk loci for melanoma 
susceptibility are the genes CDKN2A, located on chromosome 
9p21, and CDK4 on 12ql4. The CDKN2A locus encodes for both 
proteins plG''^^** and pi 4^^^ and accounts for susceptibility in 
25-40% of melanoma families [8,9]. Mutations in CDK4 are rare 
and worldwide only three families have been reported to carry 
mutations on this gene worldwide [10,11]. However, familial 
melanoma only comprises approximately 10% of all MM cases, so 
it seems likely that there are other low-penetrance polymorphisms 
with small effects but very common in the population also 
associated with the susceptibilit)' to develop MM. The MCIR gene 
is the candidate locus par excellence associated with the appearance 
of different pigmentary phenotypes as well as to the ability to 
modulate the susceptibility to develop sporadic MM [12,13]. 
However, other these low-penetrance alleles have been reported 
which belong to biological processes such as pigmentation (ASIP, 
0CA2, SLC45A2, TYR, TYRPl), immune response [IL-lfi, IL-10, 
INF-y, TNF-a), DNA repair iERCCl, ERCC2, MGMT, TERTl, 
TRFl, XRCCl, XRCC3), and metabolism (GSTMl, GSTPl, 
GSTTl), including the vitamin D receptor [12,14—19]. 

Despite the fact that investigating the molecular alterations 
involved in the pathogenesis of MM is a topic of active research, 
the advances achieved so far are still insufficient to establish a set of 
biomarkers and molecular targets which will facilitate an early 
diagnosis, predicting the risk for metastasis in melanoma patients, 
and developing more efficient therapies against this neoplasm. In 
this context, we have conducted a case-control association study 
based on 384 SNPs distributed in 88 candidate genes, of which 
314 were successfully genot\ped. These SNPs were selected by the 
combination of the information obtainc-d from proteomic analysis 
by two-dimensional electrophoresis and mass spectrometry, and 
mRNA expression arrays performed previously by our group. The 
selected SNPs were genotyped in a case-control study using three 
series of Spanish population samples (1,883 MM patients and 
1,358 disease-free controls in total). Our results showed that the 
SNPs rs6854854, located in ANXA5 and, rs6431588 in ILKAP are 
associated with MM. 

Results 

The candidate genes to be genotyped in the present work were 
selected based on proteins and gc'n(;s differentially expressed 
between melanoma cell lines and primary melanocytes. Therefore, 
previous proteomic analyses were performed on six melanoma cell 
lines (A375, Hs294T, HT-144, 1205Lu, WM793B, JSG) and four 
primary melanocytes (HEMn-LP, HEMn-MP, HEMn-DP, 
HEMa-LP) using two-dimensional electrophoresis (2D-PAGE), 
while Affymetrix Human U133A GeneChip arrays were used to 
analyze and compare mRNA expression profiles on four different 
melanoma cell lines isolated from patients' biopsies in our 
laboratory (which were named JEM, JSGl, JSG2, MJOI) and 
three benign nevi from patients' skin lesions isolated also in our 
laboratory (named FDR, JPA, RRR) (unpubhshed data). Finally, 
we selected 88 genes involved in cell growth, cell cycle and 



apoptosis, cell signaling, transcription and stress response (Table 1). 
A total of 384 SNPs were selected on the basis of their Linkage 
Disequilibrium (LD) profiles using Haploview [20] . 

In discovery Phase I, a total of 464 patients with MM and 400 
volunteer cancer-free controls were genotyped. From the initial list 
of 384 SNPs, 70 SNPs were discarded in Phase I for the following 
reasons: 17 SNPs because they could be genotyped in less than 
85% of the samples; 5 SNPs due to a low quality genotyping; 
18 SNPs were found to be monomorphic in our control 
population and patients; and 30 SNPs were out of Hardy- 
Weinberg equilibrium (HWE) after a conservative Bonferroni 
correction for multiple testing. The list of removed SNPs is 
provided in Table SI in File SI. 

Comparison of Allele Frequencies 

Therefore, the study continued with 314 SNPs (Table S2 in File 
SI), whose allele frequencies were estimated based on our control 
samples and compared to those of the HapMap CEU population 
(people with Northern and Western European ancestry). Frequen- 
cies were very similar and a high positive correlation can be 
observed in Figure 1 (R^ = 0.92), suggesting that our study 
population is sufficiently homogeneous to conduct genetic 
association studies with minor risk of population stratification. 
However, one-sample i-test gave evidence that the corresponding 
Spanish minor allele frequency (MAF) differed from that published 
in HapMap (^-value <0.05) for 8 of the 314 SNPs considered 
(rs2069502 in CDK4, rs3731239 in CDKNIA, rs2303942 in FASTK, 
rs2497 in GDI2, rs2088702 in PEBPl, rs228275 in PSAdBS, 
rsl 7800727 in RBL2, rs6586542 in RCC2) (shown as blue dots in 
Figure 1). 

Associations with MM Risk 

After a Fisher's exact test to compare allele counts between cases 
and controls, 38 SNPs located in 31 genes were associated with 
MM in the Spanish population considering a p-mlue threshold of 
0.05. Representation of -loglO p-values for the comparison of 
minor allele frequency (MAF) between the 464 MM cases and the 
400 controls are detailed in Figure 2. Detailed information on 
SNP, gene, chromosome location, MAF, odds ratio (OR), 95% 
confidence interval (95% CI) and p-value for these 37 SNPs are 
presented in Table 2. If a more restrictive p-value threshold of 0.01 
is established, 9 SNPs remain as candidates associated with MM in 
our Spanish population. And among them, only rs6854854, 
located in the intron 2 oi ANXA5, showed an association that was 
statistically significant (p-iia/ac = 4 x 1 0 '^), after applying the 
Bonferroni correction for multiple comparisons. This means that 
individuals who possess at least one AJVXA5 rs6854854 C allele are 
protected against developing MM relative to those with the 
reference genotype (OR = 0.541; 95% CI = 0.371-0.791; p- 
value= 0.0013). The Aj\XA5 rs6854854 locus was significant under 
the additive {p-value = 9.06 x 1 0~'^) and recessive genetic models (p- 
m/mc ==5.97x10"-'). 

From tiiese 38 SNPs with a p-value <0.05, we selected 15 SNPs 
with a MAF greater than 0.05 in our control samples for 
validation in Phase II, which consisted on an independent set of 
507 MM cases and 383 controls. None of the SNPs presented a 
statistically significant association with MM at this stage (Table 3). 
However, four of them had an overall p-value <0.05 when phases I 
and II were considered together (971 MM cases and 783 controls). 
These four SNPs were: rsl 3 167522 in ARC (p-aate = 0.0035); 
rs4874163 in EEFID (p-vake = 0.007 7); rs6431588 in ILKAP (p- 
= 0.0086); and rs7212835 in PSMDll (p-Mte = 0.038). As 
association studies require large numbers of samples to detect 
weak association and increasing study size typically has a large 
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Table 1. Eighty-eight candidate genes selected for SNPs genotyping. 



Category Genes 

Cell growth, cell ANXA5, ARC, AR, AURKS, AXL, SAX, BCL2U1, CDK2, CDK4, CDKNIA, 

cycle and apoptosis CTSD, ENOl, FASTK, GAS6, GSTPl, HDGF, IGFBP5, ILKAP, IMPDH2, 

MAGEDl, MAPREl, NDN, NME2, NPMl, NRAS, PEA15, PIMl, PSMB3, 



PSMD11, RAD50, RAP IS, R8L2, RCC2, RELA, SPRR2G, TGFSl, WEE! 


Cell signaling 


AR, BAG2 CALM3, CCT7, CLICl, CSNKIGZ FSCNl, GDI2, GRSZ ITGA5, 

ITGAM, MAPKAPK3, MCLl, MYD88, PDPKl, PESPl, PKM2, RACl, 

RCNl, RPSA, SIAH2, SNAIl, THBSl, TPM2 TTCl, UBE2L6, WARS, WNT5A, YWHAZ 


Transcription 


ARID5A, CRESSP, CTSPl, CTNNSl, EEFID, EOS, HDAC5, 
HTATIP JUN, PIR, RUVBLl, RUVBL2, SPl, TFDPl, PIAS3 


Stress response 


DUSPl, GPXl, PARK7, PRDXl, PRDX3, S0D2, TXNLl 



doi:1 0.1 371 /journal.pone.0095522.t001 



effect on power, these four SNPs were genotyped in 9 1 2 new MM 
cases and 581 controls in validation Phase III. 

After genotyping 1,883 patients with melanoma and 1,358 
control subjects, one SNP in intron 3 of /LK4Pgene (rs6431588) 
was associated with increased risk of developing MM in the 
Spanish population (Fisher's test p-value = 5xl0~ ). The OR was 
1.29 (95% CI 1.12-1.48; p-«fl/Me = 4xlO"'^ (Figure 3). Therefore, 
the ILKAP rs6431588 locus was significant under the additive 
genetic model (p-Mfae = 5.8 x 10~*, using the Cochrane- Armitage 
trend test). For the other three remaining SNPs p-values were not 
statistically significant. 

ANXA5 and ILKAP Expression in Human Melanoma Cell 
Lines and Primary Melanocytes 

In order to confirm an alteration in the expression levels of 
ANXA5 and ILKAP in MM, quantitative measurement of these 
genes' expression was investigated in 10 melanoma cell lines and 3 
primary melanocytes using quantitative real time PCR (Figure 4). 



Q. 
(0 




-0.2 00 02 04 0.6 

MAF HapMap_CEU samples 

Figure 1. Comparison of minor allele frequencies (MAF), 
Spanish vs. HapMap European data. The sntall and big red circles 
show the range within one and two standard deviation (SD) of the 
mean, respectively. Blue dots represent values that significantly differ 
from HapMap European data. 
doi:1 0.1 371/journal.pone.0095522.g001 



The results showed a weak reduction in the ANXA5 gene 
expression and a significant decreased gene expression of ILKAP 
(p-value <0.05 after applying Student's Z-test) in the melanoma cell 
lines studied in comparison with primary melanocytes. Although 
some individual heterogeneity was also observed as not all the 
melanoma cell lines showed a reduction in AJ\fXA5 or ILKAP gene 
expression, the trend observed coincided with that obtained 
previously using expression microarrays, in which only ILKAP, but 
not ANXA5, showed significant differences at mRNA transcrip- 
tional levels between melanocytes and the tumor cell lines studied 
(unpubhshed data). Additionally, the results obtained from 
previous 2D-PAGE assays made in our laboratory suggested a 
decrease of both proteins ILKAP and Annexin A5 in melanoma 
cell lines compared to primary melanocytes. A Western blot 
confirmed a significant decrease of Annexin A5 protein amount in 
melanoma cell lines (Figure 5). These results suggest a deregulation 
of ILKAP at the transcriptional level, whereas the deregulation 
seems to be at posttranslational level for Annexin A5. 

ILKAP Coding and Promoter Region Sequence Analysis 

As a change in the nucleotide sequence within the gene coding 
region can alter the expression and/or mRNA stability and thus 
the final protein concentration, we sequenced a cDNA fragment of 
1,259 bp covering the whole ILKAP exons (Chr 2:239,079,043- 
239,112,324 according to UCSC Genome Browser, GRCh37/ 
hgl9), in order to assess if there were nucleotide changes in the 
melanoma cell lines. However, the chromatograms showed no 
changes in the nucleotide sequence of the ILKAP. 

In the absence of coding region mutations that could be 
associated to the ILKAP expression deregulation observed in 
melanoma cell lines, we hypothesized that there could exist 
mutations in the promoter region that could be altering the 
transcriptional activity of ILKAP. Therefore, the ILKAP gene 
promoter region was amplified from genomic DNA using two 
pairs of primers that resulted in one fragment of 1,443 bp and 
another overlapping fragment of 1,110 bp. The sequences 
obtained for each sample were aligned and assembled, enabling 
the reading of a sequence of 2,069 bp (Chr 2:239,112,648- 
239,1 14,717). Twelve single nucleotide variants were found, all of 
them corresponding to already polymorphisms described. Table 4 
shows the SNPs identified in the promoter region of ILKAP and the 
genotypes that each cell line presents. The 12 polymorphisms 
observed in the tested cell lines clustered into three different 
haplotypes, which are sorted from the position Chr 2:239,1 14,642 
to the position Chr 2:239,113,023. The melanoma lines A375, 
HT-144 and JSG are homozygous for Haplotype 1 (TACCG- 
GATCCGA). The ceU lines 1205Lu, WM793B and HEMn-MP 
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Figure 2. SNP association results. The -loglO of the allelic p-values from 314 SNPs comparing 464 melanoma patients and 400 controls of 
Spanish origin at Phase I are represented. The chromosomal SNP distribution is shown. The SNP rs6854854 in ANXA5 remained statistically significant 
after Bonferroni correction (.p-value <0.00015). 
doi:1 0.1 371/journal.pone.0095522.g002 



are homozygous for Haplotype 2 (CGTGGCGCCAAG). Finally, 
the cell lines Hs294T, HEMn-LP and HEMn-DP are homozygous 
for Haplotype 3 (CGTGACGCTAAG). 

Tajima's D Test Analysis to Identify Putative Signatures of 
Positive Selection 

Due to the high frequency of polymorphisms observed in the 
promoter region ot ILKAP, we decided to use the Tajima's D test to 
investigate if this locus was subjected to the elfect of Natural 
Selection in order to add evolutionar)' relevance to the function of 
this locus. Thus, we used data from 1000 Genomes Project to 
analyze a fmal region of 66 Kb (chr2:239066043-239 132324 
according to UCSC Genome Browser website, GRCh37/hgl9), 
which intruded the promoter and the coding region of ILfuiP. The 
African population (n = 246) did not show a statistically significant 
value of Tajima's D (Z)=0.7; p-M/Me = 0.16). However, the values 
of Tajima's D obtained for Europeans {D—2.54; p-value — 0. 010; 
n = 380) and for the Asian population (Z)=2.50; p-value = 0.013; 
n = 286) were statistically signilicant. These large and positive 
Tajima's D test in suggests that balancing selection could be 

playing a role on the evolutionary history of this locus. 

Discussion 

Instead of selecting among candidate genes already described in 
the literature, the present work started by searching for genes 
whose expression (both at the protein and mRNA levels) is altered 
in melanoma cell lines as compared to normal melanocytes. After 
identifying the candidate genes in this way, a set of SNPs was 
selected for each corresponding gene. The result was the selection 
of 88 loci and 384 SNPs, of which 314 fulfilled our quality criteria. 

The allele frequencies observed in our control samples were 
highly correlated to the HapMap CEU population (R^ = 0.92). 
Therefore, our experimental results are in good agreement with 
the recent report by Gayan et al. showing that the Spanish 
population is similar to Western and Northern Europeans and 
sufficiently homogeneous to conduct genetic association studies 
with minor risk of population stratification [21]. However, we also 
saw that in our control population the MAP for 8 of the 314 SNPs 
considered differed from that published in HapMap. This may be 



due to a role of thesi; SNPs (or of really linked SNPs) in adaptive 
differences to environmental conditions because the frequency of 
these 8 SNPs seems to difiFer across populations with difiFerent 
geographic location (Table S3 in FUe 81). 

From discovery Phase I, 38 SNPs showed p-values below 0.05 
and, of these, only tlu; ANXj15 rs6854854 SNP remained 
statistically significant after Bonferroni correction. The annexins 
are a super-family of closely related calcium and membrane 
binding proteins which show cell type specific expression. Twelve 
annexins, named as annexins Al— All and A13, have been 
described common to vertebrates [22,23]. This protein family has 
a wide variety of cellular functions including vesicle trafficking, cell 
division, apoptosis, calcium signahng and growth regulation 
[24,25] . Although Annexin A5 was the first annexin characterized 
for three-dimensional structure in 1990 [26], its exact physiolog- 
ical function has not been fully understood. Recent data suggest 
that the invasion capacity, a main characteristic of tumors, is at 
least in part regulated by Annexin A5 in different cancer t\^)(;s 
[27,28,29]. Thus, Wehder et al. (2009) detected a decreased 
migration activity and invasion capability of head and neck 
squamous cell carcinoma after lacking XyVXlJ [29]. Therefore, the 
reduction of Annexin A5 protein amount observed in melanoma 
cell lines may be changing motile capacity of the tumor cells. 
However, it should be aware that ANXA5 seems to have specific 
effects on distinct types of tumors [27], so experimentally 
functional analysis are needed to determine the specific role of 
Annexin A5 in MM. 

Thereby, although the results obtained suggest that allele C in 
rs6854854 has a protective role against melanoma ip-val- 
«^=4xl0 •'), with an OR of 0.541 (95% CI 0.371-0.791; p- 
»a/Me = 0.0013), after genotyping 464 MM patients and 400 
disease-free controls, the ANXA5 rs6854854 SNP did not reach 
statistical significance in validation Phase II. It is possible that the 
difficulty in rephcation of results may be due to the existence of a 
change in the minor C allele frequency in the control group of 
Phase II. The C allele frequency in HapMap for European 
population (CEU) is 0.110 which matches the frequency obtained 
for the control group in Phase I of the study, however, the control 
population tested in the Phase II showed a minor C allele 
frequency of 0.075. An analysis of allele frequencies in difierent 
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populations using data from 1000 Genomes Project supports a 
possible population structure for rs6854854 (Table 5). These data 
suggest the need for a new validation in a population where C 
minor allele frequency remains similar to that described in 
HapMap for European population. On the other hand, although 
using HWE as a screening tool removed part of the SNPs selected 
for analysis, SNP rsl7718 (out of HW equilibrium and located in 
the AMXA5 3'UTR) turned out to be significantly associated with 
MM. Fardo et al. (2009) found that true disease susceptibility loci 
subjected to various patterns of genotype miscalls can be largely 
out of HWE and, thus, be candidates for removal before 
association testing [30]. 

After genotypig a total of 1,883 MM cases and 1,358 controls, 
SNP rs643 1 588, located in ILKAP, was a new locus associated with 
a higher susceptibility to MM. Therefore, our data suggest that 
individuals who possess the T allele in the /Z£1P rs6431588 locus 
are more likely to develop MM, fitting to an additive model of 
inheritance. 

The 1 000 Genomes database shows that the frequency of the T 
allele in rs6431588 varies significantly depending on die 
geographical region considered (Table 5). The minor allele of 
rs6431588 appears more frequently in the European populations 
than in the rest of geographic regions. In fact, it seems that 
rs6431588 has a population structure where the frequency of the T 
allele increases from Africa towards the North of Europe. Some 
findings support the hypothesis that latitudinal genetic diversity 
gradients are present in humans and reflect genetic adaptations to 
different environmental pressures that have shaped the human 
genome [31,32,33]. Latitude appears to provide a good proxy for 
the selective pressures that shaped variation in our genome 
because it is correlated with different variables like mean winter 
and summer temperatures, rainfall or ultraviolet radiation 
exposure, which appears to be the predominant environmental 
risk factor for MM. 

We also have experimeiitaUy observed that the levels of ILKAP 
gene expression decrease in melanoma cells compared to normal 
melanocytes. In this regard, the ILKAP (Integrin-linked kinase- 
associated serine /threonine phosphatase 2C) plays a role in the 
regulation of diverse processes such as cell cycle progression, 
migration and ceU death, and appears to have an important role in 
oncogenic transformation [34,35,36]. Thus, it is possible that 
decreasing ILKAP expression levels favor a constitutive activation 
of ILK (Integrin-linked kinase). In its turn, this activation 
inactivates instead GSK3[5 which favors the stabilization and 
nuclear translocation of fi-catenin, which results in the subsequent 
activation of the TCF/LEF family of transcription factors that 
promote cell survival and proliferation [35]. In fact, a high 
activation of ILK is associated with poor outcome in patients with 
melanoma [37]. Likewise, low levels of ILKAP may reduce 
apoptosis induced by Tumor Necrosis Factor alpha (TNFa) and 
the presence of reactive oxygen species (ROS), as well as reduce 
the formation of complexes with RSK2 (Ribosomal protein S6 
kinase-2) in the nucleus and consequently enhance the expression 
of Cyclin Dl (a RSK2 downstream substrate), which ultimately 
promotes tumor cell survival and proliferation [36]. Researchers 
have thus far considered ILKAP a cytoplasmic protein, however, 
its location also in the nucleus opens a window to unknown 
functions of ILKAP. Anyway, ILKAP and Annexin A5 functions 
are not the main focus of the present work and it would be 
necessary to perform experimentaUy functional analysis to 
determine whether the down-regulation of these genes has actuaUy 
a role in MM. 

On the other hand, despite the emphasis put on functional 
analyses of coding SNPs, many SNPs are located in non-coding 
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Figure 3. Forest plot showing the odds ratios and 95% confident intervals for the four SNPs most associated with melanoma risk 
predisposition. Dots represent odds ratios for each phase, discovery (Phase I) and validations (Phases II and III). The non-shown p-values are not 
statistically significant. 
doi:1 0.1 371/journal.pone.0095522.g003 



regulatory regions whose exact functions are not yet clear, but tliat 
could be influencing the binding affinity of transcription factors, 
and thus they could be exerting an important biological regulatory 
role [38]. As we detected no variation within the coding region, we 
decided to investigate the 5' region of ILKAP, where we detected 
12 SNPs. The minor allele frequencies for the 12 SNPs found in 
the ILKAP gene promoter region are very similar for all the three 
compared populations, except for rs 11 694064 and rs 11695186 
whose minor allele frequencies are similar in the European and 
African populations but are almost zero in the Asian population, 
according to the data collected from 1000 Genomes Project 
(Table 5). Although we detected three different haplotypes in the 
melanoma cell lines and primary melanocytes studied, a total of six 
different haplotypes have been reported worldwide in the 1000 
Genomes database, whose genealogical relationship suggests the 
existence of two different main lineages (Figure 6). This fact 
typically occurs under certain non-neutral conditions such as 
under balancing selection. Tajima's D test suggests that this could 
actually be the case, which adds evolutionary relevance to the 
diversity patterns of ILKAP. We ignore which functional mecha- 
nism is shaping ILKAP diversity but it is likely relevant for the 



survival of the species and adds meaning to the association of 
rs6431588 to melanoma risk. 

In summary, we have found that ANXA5 and ILKAP expression 
are down-regulated at the transcriptional level in MM cells 
compared to melanocytes, suggesting that these two genes could 
have an important role in MM. Moreover, we have detected two 
SNPs associated with MM in these genes for the first time: 
rs6854854 on the ANXA5 gene and rs643 1 588 on the ILKAP gene. 
Both SNPs show different allele frequencies among populations 
that differ in geographical location and additionally ILKAP region 
is under balancing selection, which suggest the role of the 
environment in MM susceptibility. 

Materials and Methods 

SNPs Genotyping 

Ethics statement. AH subjects gave written informed consent 
and the study was approved by the Ethics Committee of Cruces 
and Basurto Universitary Hospitals (Bizkaia, Spain); Gregorio 
Maranon Hospital (Madrid, Spain) and University Clinic Hospital 
(V alencia,Spain). 
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Genotyping study participants. A total number of 1,88,3 
patients with melanoma and 1,358 cancer-free controls were 
genotyped in three phases. Discovery Phase I: A total of 464 
patients with malignant melanoma (MM) were recruited from 



March 2004 to March 2010 at the Dermatology Ser\dces of two 
hospitals from the Basque Country: 246 from Basurto University 
Hospital (Bilbao, Spain) and 218 from Cruces University Hospital 
(Barakaldo, Spain). Similarly, 400 cancer-free controls from the 
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Basque Country were recruited. The melanoma patients and 
cancer-free controls in this study were of Cauca.sian origin based 
on their self-declared ethnicity and the overall demographics of the 
region. Written informed consents were obtained from all the 
study participants and the study was approved by both Hospital 
Ethics Committees. Validation Phase lit The second phase of 
the study, consisting of an independent validation series, was 
composed of 507 patients with melanoma recruited from the 
Dermatology Services of three Spanish hospitals: 211 from 
Gregorio Maranon General University Hospital (Madrid, Spain), 
188 from La Paz University Hospital (Madrid, Spain) and 
108 from Ramon y Cajal University Hospital (Madrid, Spain). A 
total of 383 cancer-free controls were recruited from the same 
geographical region (Madrid, Spain). Validation Phase III: The 
third phase of the study was composed of 912 patients with 
melanoma recruited from the Dermatology Services of five 
different hospitals: 92 from Basurto University Hospital (Bilbao, 
Spain), 122 from Castellon Province Hospital (Castellon, Spain), 
207 from Hospital Dr. Negrin from Las Palmas (Gran Canaria, 
Spain), 166 from Gregorio Maraiion General University Hospital 
(Madrid, Spain), and 373 from Instituto Valenciano de Oncologia 
(Valencia, Spain). Similarly, 581 cancer-free controls were 
recruited from the geographical regions covered by the hospitals 
involved in this third phase of the study. 

SNP selection and genotyping. Proteomic analyses were 
performed on six melanoma cell lines (A375, Hs294T, HT-144, 
1205Lu, WM793B, JSG) and four primary melanocytes (HEMn- 
LP, HEMn-MP, HEMn-DP, HEMa-LP). On the other hand. 



AfFymetrix Human U133A GeneChip arrays were used to analyze 
and compare niRNA expression profiles on four melanoma cell 
lines isolated from patients' biopsies in our laboratory (which were 
named JEM, JSGl, JSG2, MJOI) and three benign nevi from 
patients' skin lesions isolated also in our laboratory (named FDR, 
JPA, RRR). 

By combining the information obtained from the results of 
protein analysis by two-dimensional gel electrophoresis (2-D 
PAGE), Western blot and mRNA expression microarrays 
performed previously by our group (data no shown), we produced 
a list of candidate genes related to melanoma and selected the 
potential SNPs using the program Haploview (www.broadinstitute. 
org/haploview) and HapMap Phase 1 & 2 full dataset (http:// 
hapmap.ncbi.nlm.nih.gov/). Finally, a total of 384 SNPs located 
in 88 different genes were chosen related with cancer for their 
involvement in cell growth, cell cycle and apoptosis, cell signaling, 
transcription and stress response. 

In Phase I, SNPs were genotyped using the GoldenGate 
Genot)'ping Assay system according to the manufacturer's 
protocol (lUumina, San Diego, CA, USA) using services from 
Progenika Biopharma, Bizkaia, Spain. Genotyping was carried out 
using 350 ng of DNA per reaction and genotypes were called 
using the proprietary software supplied by lUumina (BeadStudio, 
V3.1.3.). In validation phases, 15 SNPs (rs6854854, AAX4,5; 
rsl3167522 and rsl882619, APC; rs4135385, CTNNBl; 
rs4874163, EEFID; rs7997328, GAS2; rs6431588, ILKAF, 
rs7212835, PSMDll; rsl7800727, RBL2; rsl 1719889, RUVBLl; 
rs2241715, TGFBl; rs3733868 and rs6556466, TTCl; rs655539. 
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Figure 6. Haplotypes in the promoter region of ILKAP gene in different populations obtained from 1000 Genomes Project. Although 
we detected three different haplotypes among the melanoma cell lines and the primary melanocytes studied (underlined haplotypes), a total of six 
different haplotypes have been reported worldwide. Each haplotype is represented with a circle, whose size is proportional to their frequency in the 
global population. The genealogic relationship suggests the existence of two different main lineages, showed in green and red circles. 
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TXNLl; and rsl 1042431, WEE]) were genotyped using the 
KASPAR SNP Genotyping System (KBiosciences, Hoddesdon, 
UK). The PGR was carried out according to the manufacturer's 
instructions. The genotype of each sample was determined by 
measuring aUele-specific final fluorescence in an ABI Prism 
7900HT Detection System, using the SDS 2.3 software for allele 
discrimination (Applied Biosystems, Foster Gity, CA, USA). As a 
quality control measure, one sample duplicate and a non-template 
sample per 96-well plate were included. 

Genotyping statistical analysis. For all polymorphisms 
studied in Phase I, Fisher's exact test was used to account for 
diflerences in allele frequencies between HapMap CEU data and 
population data from the Basque Gountry, to test for deviations 
from Hardy-Weinberg equilibrium (HWE) among controls, and to 
compare allele counts between cases and controls. Correction for 
multiple testing was carried out using the Bonferroni method 
based on a final set of 314 accepted marker loci. Genotype-related 
odds ratios, their corresponding 95% confidence intervals and 
associated p-values were estimated via logistic regression using 
SPSS V. 17 and the online software "Hardy-Weinberg equilibri- 
um" (http:/ /ihg2. helmholtz-muenchen.de/ cgi-bin/hw/hwal .pi). 
Statistical analysis of the SNPs located in chromosome X was 
done with the genotypes of women only. Genotype frequencies 
were also compared between MM cases and controls using the 
Cochran-Armitage trend test. 

ILKAP Expression and Sequencing in Melanomas and 
Primary Melanocytes 

Cell Unes. In the present work, ten melanoma cell lines and 
three primary melanocytes were used. The three primary human 
melanocytes were purchased from Invitrogen (Cat. No. G-002-5G 
for lightly pigmented neonatal foreskin, HEMn-LP; Gat. No. C- 



102-5G for moderately pigmented neonatal foreskin, HEMn-MP; 
and Cat. No. C-202-5G for darkly pigmented adult foreskin, 
HEMn-DP). All primary human melanocytes were grown in 
Cascade Medium 254 supplemented with Cascade Human 
Melanocyte Growth Supplement (both from Invitrogen; Carlsbad, 
CA, USA) in the absence of antibiotics. 

Likewise, we cultivated ten different melanoma cell lines. Of 
these, A375 (ATCC GRL-1619), Hs294T (ATCC HTB-140), HT- 
144 (ATCC HTB-63), WM793B (ATCC CRL-2806), and 1205Lu 
(ATCC CRL-2812) were purchased from American Type Culture 
CoUection (Rockville, MD, USA); RPMI7951 (ACG76), COLO- 
800 (AGC193), MEL-HO (AGC62), and MEL-Juso (ACG74) were 
obtained from Innoprot (Derio, Bizkaia, Spain); and JSC was 
established and characterized in our laboratory from a surgical 
primary melanoma as described previously [39]. The melanoma 
cell lines were cultured in appropriate medium supplemented with 
10% fetal bovine serum (FBS), 2 mM L-glutamin and antibiotics 
according to the manufacturer's instruction. All primary human 
melanocytes and melanoma cell lines were cultured at 37°C with 
5% CO2 and 95% humidity. 

Gene expression by quantitative real-time PGR (RT- 
qPCR). Total RNA was isolated from ten melanoma cell lines 
(A375, Hs294T, HT-144, 1205Lu, WM793B, JSG, MEL-HO, 
MEL-Juso, COLO-800, RPMI7951) and three primary melano- 
cytes (HEMn-LP, HEMn-MP, HEMn-DP) using the RNeasy Mini 
kit (Qiagen Inc., HUden, Germany). For each sample, cDNA was 
synthesized from 1 |j,g total RNA using the iScript''" cDNA 
Synthesis kit (Bio-Rad Laboratories, Hercules, GA, USA) accord- 
ing to the manufacturer's instructions. Real-time RT-PGR assays 
were carried out using an iCycler PGR platform (Bio-Rad 
Laboratories, Hercules, GA, USA). The reaction mixture 
contained 0. 1 (il cDNA from the reverse transcription reaction. 
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together with forward and reverse specific primers and iQ_ 
SYBR® Green Supermix (Bio-Rad Laboratories, Hercules, CA, 
USA) in a final reaction volume of 20 ^1. The PGR reaction began 

by heating at 95°C for lOmin, followed by 45 cycles of 
denaturation at 95 °G for 30s, annealing at the corresponding 
temperature for each gene (56-6 1°G) for 20s and extension at 
72°G for 30s. Each assay included a negative control with no 
template. Expression data were generated from 2 amplification 
reactions with samples and controls run in triplicate. Optical data 
obtained by real-time PGR were analyzed using the MyiQ, Single- 
Golor Real-Time PGR Detection System Software, Version 1.0 
(Bio-Rad Laboratories, Hercules, CA, USA). Melt Gurve analysis 
of each PGR assay and 1.5% agarose gel electrophoresis analysis 
of randomly selected samples were performed to confirm the 
specificity of the amplification products. To normalize expression 
data obtained from the studied genes, we used the expression of 
three different housekeeping genes {ACTB, GAPDH, and RPS15) 
and the G(;ne Expression Macro Software v. 1.1 (Bio-Rad 
Laboratories, Hercules, GA, USA), where the relative expression 
values were computed by the comparative Gt method [40,41]. 
The sequences of primers used were: AMXA5, forward 5'- 
GAGGGGATGTTGGTGGTTG-3' and reverse 5'-GAGGGT- 
GAAATAAAGCCTGAG-3'; ILKAP, forward 5'-AAGTTTG- 
TAAAGGGTGTTGGGTG-3' and reverse 5'-GTGGGTGAT- 
GTGGTTCAGGAG-3'; ACTB, forward 5'-AGATGAGGCA- 
GATGATGTTTGAG-3' and reverse 5'-GTGAGGGGAGTC- 
GATGAGG-3'; GADPH, forward 5'-GGTGTTGGAGAGTGA- 
GGCG-3' and reverse 5'-CGACGAAATGCGTTGACTCG-3'; 
RPS15, forward 5'-TTGGGCAAGTTGACCTAGG-3' and re- 
verse 5'-GGGGGGGGGGATGCTTTACG-3'. 

ILKAP sequencing. ILKAP was sequenced in ten melanoma 
ceU lines (A375, Hs294T, HT-144, 1205Lu, WM793B, JSG, 
MEL-HO, MEL-Juso, COLO-800, RPMI7951) and three 
primary melanocytes (HEMn-LP, HEMn-MP, HEMn-DP). 
DNA was amplified by PGR with specific primers and using th(- 
reaction mix ImmoMix^^ Red (Gentaur, Kampenhout, Belgium) 
according to the following protocol: a lOmin denaturation at 
95°G, 35 three-step cycles (95°G for 30s, 56-60°G for 30s, and 
72°G for Imin), and lOmin at 72°G in an iGycler PGR platform 
(Bio-Rad Laboratories, Hercules, GA, USA). The removal of the 
unincorporated deoxynucleotide triphosphates and primers was 
performed using High Pure PGR Product Purification kit (Roche 
Molecular Biochemicals, Madrid, Spain). The purified DNA and 
3.2 pmol of either the forward or reverse primer were used in 
standard cycle sequencing reactions with an ABI PRISM BigDye 
Terminator kit and run on an ABI PRISM 310 genetic analyzer 
(both PE Applied Biosystems, Foster Gity, GA). The analysis and 
alignment of sequences were performed using Ghromas software 
and the BioEdit Sequence alignment Editor and the reference 
sequence from UGSG Genome Browser website, GRGh37/hgl9 
(http://www.genome.ucsc.edu/) and NGBI dbSNP 138 (http:// 
www.ncbi.nlm.nih.gov/ snp/). The s(;quen('es of primers used for 
ILKAP promoter region sequencing were: first pair, forward 5'- 
TGTTTGTGTGCCCATCAAGC-3' and reverse 5'- 
ATTCTGGGCAATTTCGATGA-3'; second pair, forward 5'- 
TTGGAAGGGTGGAATAAAGG-3' and reverse 5'- 
TTGTGGAGGTGTTGGGATGT-3'. The sequences of primers 
used for the sequencing of ILKAP exons were: forward 5'- 
TGAGTGTGTGTGGGTGGTG-3' and reverse 5'-AAGTGAA- 
TAGGATGGGTGG-3'. Genomic DNA was used for ILKAP 
promoter region sequencing, while cDNA was used for the 
sequencing of ILKAP exons. 

Tajima's D test. The variation of nucleotide patterns from 
the neutral expectation was tested by the Tajima's D test using 



DnaSP v.5.10.1. Gene diversity is controlled by the parameter 
theta (9 = 4 jV^ ^i, where jV^ is the efiective population size and \l 
the per generation mutation rate). Several sample- based 

estimators of theta (9) exist, all based on the site-frequency 
spectrum of the mutations (SFS), that is, the distribution of the 
proportion of sites where the mutant is at frequency x. Tajima's D 
test compares 9k and 9^ asking about the occurrence of rare and 
common variants [42] . This test takes into account the number of 
nucleotide positions at which a polymorphism is found or, 
equivalentiy, the number of segregating sites, k, and the average 
per nucleotide diversity, n. Using some mathematical expressions, 
if the nucleotide sequence variation among our haplotypes is 
neutral and the population from which we sampled is in 
equilibrium with respect to drift and mutation, then Tajima's D 
test should be indistinguishable from zero. If it is either negative or 
positive, we can infer that there's some departure from the 
assumptions of neutrality and/ or equilibrium. 

A region of 66 Kb (chr2:239,066,043-239, 132,324 according to 
UGSG Genome Browser website, GRGh37/hgl9), including the 
promoter and coding region of ILKAP gene, was analyzed and 
compared among European, Asian and African populations from 
1000 Genomes Project. Standard coalescent simulations, as 
implemented in DnaSP, were used to estimate the statistical 
significance of the D values. 

Western blot. Melanoma cells (A3 75, Hs294T, HT-144, 
1205Lu, WM793B, JSG, MEL-HO, MEL-Juso, GOLO-800, 
RPMI7951) and primary melanocytes (HEMn-LP, HEMn-MP, 
HEMn-DP) were harvested by trypsinization, washed with PBS 
and lysed in RIPA lysis buHer (80 mM Tris-HGI pH 8, 150 mM 
NaGl, 1% NP 40, 0.5% sodium deoxycholate, 0.1% SDS) 
containing Protease Inhibitor GocktaU (Sigma-Aldrich Quimica, 
S.A., Madrid, Spain) for 15 minutes on ice. Lysates were then 
cleared by centrifugation at 10,000 g for 5 minutes and total 
protein concentration was determined. Fifty micrograms of total 
proteins from each sample were resolved by electrophoresis on an 
SDS-polyacrylamide gel and then transferred onto a nitrocellulose 
membrane (Whatman GmbH, Dassel, Germany). The blots were 
incubated with PBS containing 5"/() non-fat milk and 0.1% Tween- 
20 for 1 hour to block nonspecific binding, and then incubated 
with an appropriate dilution of primary antibody at 4°G for 
overnight. The primary' antibodies used were anti-human Annexin 
A5 (ab54775) and y-tubulin (abll320) antibodies (Abeam, Inc, 
Gambridge, GA). After washing, membranes were incubated for 
1 hour with horseradish peroxidase-hnked secondary antibody. 
Finally, proteins were visuafized by enhanced chemUuminescence 
using the SuperSignal® West Pico GhemHuminescent Substrate 
(Thermo Scientific, Rockford, IL, USA) and the intensity of each 
band was measured using ImageJ software. 

Supporting Information 

File SI Table SI. SNPs removed in the discovery Phase I of the 
genotyping study. Table S2. List of 314 successfully genotyped 
SNPs, HapMap_GEU MAF, Spanish MAF, and HWE p-»ate. 
Table S3. Minor allele frequency in different populations for the 
8 SNPs appearing as outiiers in Figure 1 (Data from 1000 
Genomes Phase I May 201 1). 
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